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Abstract 

A  general  scheme  for  the  VLSI  implementation  of  auditory  wavelet  transforms  is  proposed 
using  switched-capacitor  (SC)  circuits.  SC  circuits  are  well  suited  for  this  application  since 
the  dilation  constant  across  different  scales  of  the  transform  can  be  precisely  implemented  and 
controlled  by  both  capacitor  ratios  and  the  clock  frequency.  The  hardware  implementations  are 
made  possible  by  several  new  circuit  designs.  Specifically,  extremely  area-efficient  designs  are 
presented  to  implement  very  large  time-constant  filters  such  as  those  used  to  process  speech  and 
other  acoustic  signals.  The  designs  employ  a  new  charge  differencing  technique  to  reduce  signif¬ 
icantly  the  capacitance  spread  ratios  needed  in  the  filter  banks.  Also,  a  new  sum-gain  amplifier 
(SGA-SI)  is  designed  which  permits  several  inputs  to  be  sampled  with  the  same  phase.  The  pro¬ 
posed  circuits  have  been  fabricated  using  a  2/xm  CMOS  double-poly  process.  Preliminary  data 
and  performance  measures  of  the  circuits  are  very  encouraging  and  are  presented.  Two  possible 
architectures  for  implementing  the  wavelet  transform  are  discussed  and  compared:  parallel  and 
cascade  filter  banks.  Responses  of  both  filter  banks  are  simulated  using  SWITCAP-II.  Finally, 
we  shall  also  briefly  discuss  the  utility,  from  an  implementation  point  of  view,  of  decomposing 
the  transfer  functions  of  the  filter  banks  into  rational  form  using  a  recently-developed  wavelet 
system  (WS)  technique. 


1  Introduction 

Wavelet  transformations  have  found  a  wide  range  of  applications  in  the  processing  of  signals  such 
as  speech  and  images  [1],  In  biological  systems  such  as  the  auditory  and  visual  nervous  systems. 
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affine  wavelet  transforms  serve  as  excellent  models  of  the  operations  leading  to  the  generation  of 
robust  and  perceptually  accurate  representations  of  the  input  signals[2].  However,  a  principal  ob¬ 
stacle  to  the  wider  utilization  of  these  and  other  nonorthogonal  wavelet  based  algorithms  is  the 
heavy  computational  cost  of  the  wavelet  transform  stage.  Consequently,  hardware  implementations 
have  been  an  attractive  option  to  achieve  real-time  performance.  For  instance,  a  number  of  at¬ 
tempts  at  hardware  implementations  of  the  cochlear  filter  bank  have  employed  analog  designs  using 
subthreshold  operational  transconductance  amplifiers  to  construct  a  cascade  of  second-order  filter 
stages  [3,  4,  5].  Although  successful  in  many  respects,  these  designs  suffer  from  the  effects  of  par¬ 
asitic  capacitances  and  nonuniformity  of  the  fabrication  process,  and  thus  require  post-fabrication 
tuning  to  compensate  for  these  effects.  To  alleviate  some  of  the  difficulties  associated  with  such 
analog  designs,  we  consider  a  different  approach  which  employs  a  bank  of  switched-capacitor  filters 
(SCF’s)  [6,  7].  SCF’s  in  general  have  extremely  precise  and  reliable  response  characteristics  that 
would  obviate  the  need  for  any  post-fabrication  tuning. 

A  key  feature  of  using  SC  circuits  for  implementing  wavelet  transforms  is  that  dilations  of  a  given 
filter  may  be  easily  and  very  precisely  controlled.  Two  distinct  mechanisms  exist  for  implementing 
dilations  via  SC  circuits.  Given  a  filter  G(s),  and  its  dilated  version  G(as),  the  direct  approach  of 
SC  circuit  implementation  permits  control  of  the  dilation  constant  a,  to  within  0.1%  error,  since  a 
depends  only  on  the  ratios  of  certain  capacitors  and  not  their  absolute  values.  A  second,  and  perhaps 
more  interesting  method  of  controlling  a,  in  a  SC  circuit  implementation  involves  controlling  the 
various  clock  frequencies  of  the  circuits.  This  second  approach  implementing  dilations  is  discussed 
further  in  Section  5.  Such  accuracy  is  in  general  unachievable  using  conventional  analog  designs. 
The  circuits  described  here  are  very  area-efficient,  compensated  for  non-ideal  effects,  and  are  free 
from  the  effects  of  parasitic  capacitances. 

Several  difficulties  arise  in  designing  switched-capacitor  filters  to  implement  wavelet  transforms 
of  acoustic  signals.  The  most  serious  is  the  need  for  a  frequency  range  which  is  broad  and  stretches 
to  relatively  low  frequencies.  Conventional  switched-capacitor  circuit  designs[8,  9]  require  a  capac¬ 
itance  spread  ratio  of  approximately  l/(f IqT),  where  fio  is  the  pole  frequency  of  the  linear  section 
or  the  second  order  filter  (also  known  as  a  biquad),  and  T  is  the  sampling  period.  Thus,  for  a  filter 
processing  low  acoustic  frequencies,  this  ratio  becomes  very  large,  hence  requiring  a  large  silicon 
area.  Another  problem  with  existing  SCF’s  is  that  area  efficient  designs  are  achieved  by  signal 
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attenuation,  rendering  them  vulnerable  to  op-amp’s  non-ideal  effects  such  as  input  offset  voltage 
and  finite  DC  gain[10]. 

In  order  to  overcome  these  difficulties,  new  system  architectures  and  SC  circuit  designs  are 
presented  that  would  facilitate  the  VLSI  implementation  of  arbitrary  wavelet  transforms  as  parallel 
or  cascade  filter  banks.  In  order  to  illustrate  the  design  concepts,  a  specific  implementation  of  a 
32-channel  wavelet  transform  is  carried  out.  The  filter  bank  mimics  the  transformations  observed 
in  the  cochlea  of  the  inner  ear  which  are  discussed  in  detail  in  [11].  In  the  following  section,  we 
discuss  two  overall  system  design  considerations:  (1)  the  specification  of  the  filter  shapes,  and  (2) 
the  architecture  of  the  system,  specifically,  the  advantages  and  disadvantages  of  parallel  versus 
cascade  implementations  of  the  filter  bank.  In  Section  3,  we  discuss  in  detail  the  design  of  the 
new  SC  circuits  needed  to  implement  the  parallel  filter  bank.  Examples  of  such  circuits  are  a  new 
very  large  time-constant  (VLT)  switched-capacitor  circuit  employing  gain  and  offset  compensation 
(GOC)  and  charge  differencing  (CD)  techniques  [12],  and  a  sum-gain  amplifier  (SGA-SI)  which 
permits  inputs  to  be  sampled  with  the  same  phase  to  eliminate  the  need  for  sample-and-hold 
circuitry  at  the  output  of  each  individual  filter  (or  filter  section).  The  performance  of  IC  fabricated 
versions  of  these  circuits  is  also  discussed  briefly.  In  Section  4,  the  cascade  filter  bank  option  is 
considered.  Finally,  simulation  results  of  the  two  filter  bank  architectures  are  also  compared. 

2  Overall  System  Design  Considerations 

2.1  Design  of  Filter  Transfer  Functions 

In  the  design  of  a  wavelet  transform,  the  first  issue  to  be  tackled  is  that  the  exact  shape  and  number 
of  filters  to  be  used.  This  in  general  involves  application  dependent  considerations.  For  example, 
the  desired  frequency  response  of  a  single  channel  of  the  cochlear  filter  bank  may  be  empirically 
obtained,  analytically  derived  from  mathematical  models  of  the  basilar  membrane,  or  manually 
designed  to  suit  a  particular  application.  Often  the  desired  frequency  responses  take  the  form 
of  unparameterized  models  (e.g.  empirically  obtained  responses)  or  nonrational  transfer  function 
models  (e.g.  responses  derived  from  continuum  mechanical  models  of  the  basilar  membrane).  For 
the  purpose  of  circuit  implementations,  it  is  often  convenient  (as  in  the  case  of  the  parallel  filter 
bank)  to  look  for  finite-dimensional  linear  filters  to  approximate  the  desired  responses  of  the  cochlear 


filters.  It  is  well-known  that  a  linear  filter  has  a  finite-dimensional  state  space  realization  if  and 
only  if  the  corresponding  transfer  function  is  rational.  Hence  the  problem  is  one  of  constructing 
rational  approximants  to  the  desired  cochlear  filter  transfer  functions. 

Systematic  methods  of  rational  approximation  have  received  considerable  attention  in  the  con¬ 
trol  theory  and  system  identification  literature.  In  Appendix  A  we  discuss  the  use  of  the  recently- 
proposed  wavelet  system  (WS)  [13,  14]  technique  for  constructing  rational  approximants  to  the 
cochlear  filter  transfer  functions.  Our  use  of  this  technique  in  the  present  paper  is  motivated  by  the 
fact  that  the  WS  approximation  method  is  particularly  well-suited  to  approximate  systems  with 
transfer  functions  which  are  well-localized  in  time-frequency.  Secondly,  approximants  generated 
by  the  WS  methodology  are  suitable  for  hardware  implementation  via  the  SC  circuits  discussed  in 
Section  3.  The  basic  idea  behind  WS  approximations  is  the  decomposition  of  a  class  of  transfer 
functions  via  dilations  and  complex  translations  of  a  single  real-rational  function.  Specific  trunca¬ 
tions  of  such  decompositions  are  then  used  as  rational  approximants.  We  include  a  brief  description 
of  the  WS  methodology  in  Appendix  A.  For  a  detailed  treatment  and  comparison  with  methods 
based  on  the  classical  Laguerre  filters,  we  refer  to  [14]. 

We  should  emphasize  that  while  the  WS  formalism  provides  a  systematic  methodology  for 
rational  approximation  of  fairly  complicated  filter  shapes,  its  use  is  strictly  limited  to  the  parallel 
filter  implementations.  This  is  because  the  cascade  implementation  of  a  filter  bank  imposes  implicit 
restrictions  on  the  form  of  the  filter  transfer  functions  which  are  not  easily  derived  from  their  rational 
form.  In  the  parallel  filter  bank,  each  channel  can  in  principle  be  implemented  independently  of  all. 
others,  and  hence  a  rational  form  is  immediately  useful.  However,  in  an  especially  efficient  form  of 
the  parallel  bank  in  which  filter  elements  are  shared  across  channels  (see  Section  3),  we  shall  see 
that  in  general  the  rational  form  acts  merely  as  a  guide  for  the  design  process. 

2.2  Parallel  versus  Cascade  Architectures 

In  general,  there  are  two  possible  architectures  for  implementing  the  wavelet  transform:  as  a 
parallel  or  as  a  cascade  filter  bank.  Each  of  these  approaches  has  its  advantages  and  disadvantages, 
somewhat  depending  on  the  specific  application  at  hand.  The  primary  advantages  of  the  cascade 
filter  bank  is  the  ability  to  obtain  high  frequency  resolution  in  the  low  frequency  ranges  with 
lower  Q  circuits.  This,  combined  with  the  fact  that  no  sum-gain  amplifiers  are  required  makes 
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this  architecture  more  silicon  area  efficient  than  a  comparable  parallel  bank.  However,  a  host 
of  potential  disadvantages  may  balance  out  these  benefits.  To  start  with,  from  a  circuit  design 
viewpoint,  a  cascade  implementation  requires  more  robust  designs  since  an  accurate  DC  unity  gain 
and  low  offset  voltages  are  needed  at  each  stage  of  the  filter  bank  so  as  to  maintain  an  adequate 
system  dynamic  range.  The  other  obviously  undesirable  feature  of  the  cascade  filter  bank  is  that 
the  failure  of  one  filter  stage  will  affect  all  succeeding  filter  stages.  From  a  theoretical  point  of  view, 
the  cascade  filters  are  interdependent,  and  hence  one  may  not  be  able  to  design  one  filter  without 
affecting  many  others.  As  such,  the  ability  to  find  a  rational  form  for  the  filter  transfer  function 
in  a  given  tap  is  of  no  direct  consequence,  but  rather  is  only  useful  indirectly,  e.g.,  as  an  indicator 
of  the  needed  Q’s.  More  problematic  is  the  fact  that  in  implementing  a  wavelet  transform,  the 
dilation  relationship  among  the  different  channels  is  difficult  to  maintain,  resulting  in  inevitable 
implementation  errors  that  may  be  significant  in  some  applications. 

The  parallel  filter  bank  offers  significantly  greater  flexibility  since  its  filters  in  general  are  im¬ 
plemented  independently  of  each  other.  Thus,  a  rational  form  of  the  transfer  function  of  a  seed 
filter  can  be  implemented  directly  as  SCF  circuits,  with  other  filters  being  simple  dilations  of  it 
(e.g.,  through  SCF  clock  frequency  shifts).  It  is,  however,  also  possible  to  implement  more  efficient 
designs  by,  for  instance,  sharing  biquad  elements  between  adjacent  channels.  In  this  case,  the  bi¬ 
quad  elements  themselves  need  to  be  dilations  of  each  other,  i.e.,  it  is  essential  that  each  filter  be 
implemented  with  constant  Q  biquad  circuits.  As  in  the  cascade  case,  sharing  the  filter  elements 
creates  inter-filter  dependencies  that  render  the  rational  forms  of  the  filter  transfer  functions  only 
indirectly  useful  as  a  guide  for  more  general  curve  fitting  procedures.  We  shall  discuss  in  detail  an 
example  of  such  an  architecture  in  Section  3. 

3  Parallel  Filter  Bank  Implementations 

We  focus  in  this  section  on  the  design  of  a  32-channel  parallel  filter  bank  with  shared  elements. 
This  is  a  special  case  of  the  general  parallel  architecture  which  has  particularly  efficient  silicon 
area  utilization.  Generalizing  the  design  to  the  simpler,  but  less  efficient,  non-shared  filter  bank 
is  straightforward.  All  SCF’s  designed  to  implement  the  filter  biquads  are  illustrated  in  detail, 
together  with  simulations  of  their  response  characteristics.  The  circuits  are  generic  in  nature,  and 
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Figure  1:  System  topology  of  a  wavelet  transform  using  a  shared-component  parallel  filter  bank, 
can  be  readily  used  in  any  arbitrary  wavelet  transformation. 

3.1  System  Design:  Dilating-functions  Filter  Bank 

The  system  topology  of  the  shared  parallel  filter  bank  implementation  is  shown  in  Fig.  1.  In  this 
specific  design  example,  each  channel  is  implemented  by  connecting  3  adjacent  lowpass  biquads 
in  parallel.  The  output  signals  from  these  circuits  are  amplified  and  summed  together  through 
the  sum-gain  amplifier  (SGA)  to  approximate  the  desired  filter  transfer  function.  Note  in  general, 
the  number  of  biquads  used  in  each  filter  is  purely  dependent  on  the  order  of  the  filters  desired. 
The  biquads  are  designed  with  progressively  dilated  impulse  responses  which  permits  sharing  them 
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among  channels1.  In  the  specific  example  of  Fig.  1,  a  n-channel  filter  bank  only  requires  n+2 
biquads  and  n  sum-gain  amplifiers.  The  transfer  function  of  the  n-th  biquad  is  chosen  as 

Hn(s)  =  H(an~1s)  (1) 

where,  H(s)  is  a  2nd  order  lowpass  filter  in  the  s-domain,  a  is  a  dilation  constant,  and  n  =  1,  ■  •  34 

in  this  design  example.  Using  these  biquads,  the  transfer  function  of  the  m-th  channel  in  the  z- 
domain  is  approximated  by 

t=3 

Gm{z)  =  '52kiHi+m- l(z)  (2) 

1  =  1 

where,  fc;’s  are  gain  factors.  Since  the  same  set  of  fc,’s  is  used  in  every  channel,  only  one  sum-gain 
amplifier  has  to  be  designed.  The  rest  will  just  be  its  duplicates.  The  output  signals  from  biquads 
are  available  simultaneously  in  this  system  topology.  Therefore,  a  new  sum-gain  amplifier  (SGA- 
SI)  which  permits  the  inputs  to  be  sampled  in  the  same  phase  is  proposed  here.  The  new  design 
requires  no  sample- and- hold  circuitry  to  connect  the  outputs  of  a  parallel  bank  of  biquads  to  the 
input  of  a  SGA-SI,  thereby  saving  considerable  silicon  area.  The  pre-emphasis  highpass  filter  is 
often  useful  in  the  processing  of  acoustic  signals  such  as  speech  [15].  Since  the  filter  is  shared  among 
many  channels,  the  driving  capability  of  its  op-amp  has  to  be  considered  carefully. 

In  the  following  subsections,  we  shall  discuss  in  detail  the  SCF  circuits  that  we  designed  in 
order  to  implement  the  filter  bank.  In  the  Section  3.2,  we  focus  on  the  SCF’s  needed  in  the  low 
acoustic  frequency  ranges,  the  so-called  VLT  integrators  and  biquads.  New  designs  are  formulated 
here  in  order  to  produce  area-efficient  and  robust  filters.  For  the  midfrequency  ranges,  conventional 
biquads  can  be  used  as  in  Section  3.3.  The  highpass  pre-emphasis  filter  design  and  the  sum-gain 
amplifiers  are  presented  in  Sections  3.4  and  3.5,  respectively.  Finally,  simulations  of  the  filter  bank 
operation  using  all  these  SCF  circuits  are  shown  in  Section  3.6. 

3.2  Circuit  Design  I:  Area-Efficient  GOC  VLT  Biquad 

It  is  known  that  the  performance  of  switched-capacitor  circuits  can  be  degraded  by  the  non¬ 
ideal  effects  of  op-amps  such  as  offset  voltages  and  finite  DC  gains.  Since  these  effects  can  usu¬ 
ally  be  improved  simultaneously,  the  circuits  with  gain-  and  offset-compensation  are  called  GOC 

'Note  that  sharing  of  biquads  need  not  be  restricted  to  adjacent  channels,  and  that  more  general  shared-component 


topologies  are  also  possible. 


Figure  2:  GOC  charge  differencing  inverting  integrator. 

circuits[16,  17].  In  this  section,  we  describe  new  GOC  area-efficient  biquads  which  employ  a  tech¬ 
nique  called  charge  differencing[12].  These  circuits  are  specifically  designed  for  use  in  the  low 
frequency  channels  where  very  large  time-constants  (VLT),  and  hence  large  capacitances  and  sili¬ 
con  areas,  are  typically  needed  in  conventional  designs.  First,  we  shall  describe  a  VLT  integrator 
circuit,  and  then  use  it  to  design  the  larger  2nd  order  (biquad)  circuits. 

3.2.1  Gain  and  Offset  Compensated  (GOC)  VLT  Integrator 

The  GOC  inverting  integrator,  where  Col  >  Ca2,  is  shown  in  Fig.  2.  Its  operational  principle  can  be 
stated  as  follows  (assume  op-amp  is  ideal  first).  When  switches  1  are  closed,  the  charge  Vin{n)C\  is 
accumulated  in  the  capacitors  Ca  and  Ca2.  The  output  voltage  is  sampled  by  Ca i  simultaneously, 
i.e.,  the  charge  [I/!n(n)CaiCi]/(Cyi +  Ca2)  is  transferred  into  Cai-  When  the  switches  2  are  closed,  a 
charge  Vin(n)C\  is  effectively  pulled  back  to  ground  from  capacitors  Ca  and  Cai-  Since  Ca  received 
the  charge  \Vin(n)C  ACi\l  (C a+C  in  the  previous  phase,  Cai  has  to  compensate  the  net  difference 
charge  [Vin(n)Ca2Ci]/(CJ4  +  Ca2)  before  it  redistributes  the  charge  received  in  the  previous  phase 
with  C a-  The  transfer  function 

(  ,  Vout(z)  C,{Cal-Ca2)z-1!2 

M  Vin(z)  (CA  +  Cal)(CA  +  Ca2)(l-z~') 

can  be  obtained. 
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The  charge  differencing  technique  is  based  on  both  the  difference  of  capacitors  and  ratio  of 
capacitors,  thus  the  capacitance  spread  ratio  can  be  made  very  small.  Note  that  the  difference 
between  the  capacitors  Ca\  and  Ca 2  should  not  be  made  arbitrarily  small  because  of  sensitivity 
problems[12].  This  integrator  is  stray-insensitive,  glitch-free,  and  only  one  medium  capacitor  is 
needed. 

Now,  let  the  offset  voltage  be  Vos.  During  the  phase  2,  the  LHS  plate  of  capacitor  C0  is  connected 
to  ground,  thus  the  voltage  —  Vos  is  stored  in  CQ.  During  the  phase  1,  only  the  RHS  plate  of  C0  is 
only  connected  to  the  inverting  terminal  of  the  op-amp,  therefore  the  charge  held  in  C0  is  preserved. 
Consequently,  the  LHS  plate  of  C0  acts  like  the  virtual  ground,  and  this  can  be  used  to  improve  the 
non-ideal  effects  of  the  op-amp.  The  arrangement  of  switching  phases  surrounding  the  capacitor 
C0  is  very  important.  In  the  above  case,  the  virtual  ground  exists  in  the  phase  1.  If  the  switching 
phases  surrounding  the  capacitor  C0  are  interchanged,  the  virtual  ground  will  appear  in  the  phase 
2  and  this  alters  dramatically  the  degree  of  compensation.  Using  time  domain  analysis,  the  GOC 
integrator  output  is  given  by 


Lout(^  4"  )  — 


2 


+ 


Ci(Cal  -  Ca2) 

(Ca  +  Cal)(CA  +  Ca  2) 

Ca'2  ,  Cl(Ca2  ~  Cal 


+ 


Ut'n(U')  -(-  Vout(n  ) 


ri  Vo.- 


(4) 


C  a  +  Ca  2  ( C  A  +  Cal)(CA  +  Ca  2) 

For  the  non-GOC  case  (the  junction  of  C\  and  CA  is  connected  to  the  inverting  terminal  of  the 
op-amp  directly,  and  omit  the  capacitor  Ca  and  its  associated  switches),  it  can  also  be  shown  that 
the  effective  offset  voltage  is 


( 


Cal 

C  A  +  Cai 


+ 


Ca2 

C  a  +  Ctt2 


)vos. 


(5) 


Comparing  eqns.  4  and  5,  the  effective  offset  voltage  is  reduced  by  at  least  a  factor  of  2  in  the  GOC 
design.  In  the  special  case  of  Ca2  =  0,  the  output  offset  voltage  of  the  GOC  VLT  integrator  is 
approximately  inversely  proportional  to  CA2.  That  is,  the  effective  offset  voltage  is  further  reduced, 
and  with  lower  circuit  sensitivity[12,  18,  7,  19]. 


3.2.2  Gain  and  Offset  Compensated  VLT  Biquads 

With  the  above  GOC  charge  differencing  integrator,  two  types  of  biquads  can  be  built  (Figs.  3  and 
4).  Type-I  biquad  has  a  somewhat  simpler  design  equations  and  is  suitable  in  high  Q  applications, 
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Figure  5:  Signal  flow  graph  of  Type-I  GOC  charge  differencing  biquad. 


Figure  6:  Signal  flow  graph  of  Type-II  GOC  charge  differencing  biquad. 
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while  Type- II  is  more  appropriate  for  low  Q  filters.  From  the  signal  flow  graphs  shown  in  Figs.  5  and 
6  (assuming  ideal  op-amps  and  no  offset  voltages),  and  apply  Mason’s  rule,  the  transfer  function 
of  biquads  is  given  by 


TJ(  \  _  etc  1  -  (2  -  7e  -  rjc)z  1  +  (1  -  T]c)z  2 
[Z)  ~  C2  lc  1  -  (2  -  ac  -  ft)*-1  +  (1  -  ft)*"2  ’ 

where, 

=  C1Cs(Cal  -  Ca2)(Cbl  -  Cb2) 

7c  C6(CA  +  Cal)(CA  +  Ca2)(CB  +  Cbl ) 
CxC7{Cbl  -  Cb2) 

Vc  Ce(CA  +  Ca2)(CB  +  Cbly 


For  the  Type-I  biquad, 


Qr.  = 


C\C2(Cal  ~  Ca2){.Cb\  —  Ch2) 


ft 


(CA  +  Cal){CA  +  Ca2){CB  +  cbl\cB  +  cb2  +  C4) 
cycn  -  Cb 2) 


( cB  +  cbl)(cB  +  cb2  +  c4y 


(6) 

(7) 

(8) 

(9) 

(10) 


For  the  Type-II  biquad, 

C\C2(cai  -  ca2ycbl  -  cb2) 

ac  (CA  +  Cal)(CA  +  Ca2)(CB  +  Cbl)(CB  +  Cb2) 

CiC3(C6i  -  cb 2) 

o cA  +  ca2ycB  +  cbl)(cB  +  cb2y 

Consider  the  following  transfer  function  in  the  s-domain, 

rr  /  s  _  1  +  s/(Ql9.l)  +  S2/Q,l 

n^s)-  i  +  s/(Q0n0)  +  syni 


(11) 

(12) 


(13) 


Defining  aq  =  2/f and  performing  the  bilinear  transformation  (if  necessary,  the  prewarp  has  to 


be  done  before  taking  the  transformation),  the  z-domain  transfer  function  is  given  by 

rr,  s  _  {x\  +  Z1/Q1  +  1)  -  (2x1  -  2)z~1  +  (X1  -  xl!Q\  +  l)z~2 

2  ~  (Xq+  Xq/Qo  +  1)  -  (2a:2  -  2 )z~l  +  (xl  -  Xq! Qq  +  l)*-2  ' 


(14) 


Let 


a  = 

4 

(15) 

xo  T  (u.'o/ Qo)  +  1 

(3  = 

2x0/Qo 

(16) 

Xq  +  [xq!  Qo)  +  1 

7  = 

4 

(17) 

xt  +  (xi/Qi)  +  1 

n  = 

2  x1/Q1 

(18) 

x\  +  (xi/Qi)  +  1’ 

then  the  transfer  function  becomes 


H(Z)  =  ka  1  ~  i'2  ~  ~  I)?'1  +  i1  ~ 'l)z  2 

V  71-(2-q-/?)z-1  +  (1-/?)z-2' 

Let  the  capacitor  values  be  assigned  as: 

CX  =  C2  =  C0l  =  Co2  =  1 
Ca2  —  Cb2  —  2/^0 
Cai  —  Cb i  =  y  +  6  >  0 
CA  —  Cb  —  K  >  0 


(19) 


(20) 

(21) 

(22) 

(23) 


where  y  is  chosen  to  be  smaller  than  1  in  the  VLT  applications.  Comparing  the  coefficients  of 
eqns.  6  and  19,  the  K  for  the  Type-I  biquad  becomes  the  solution  of 

K 2  +  (2 y  +  S)K  +  (: y 2  +  yS  -  -£=)  =  0.  (24) 

For  flo T  <  1,  and  after  some  simplifications, 

K  *  \J 6/(Q.qT)  -  1  (25) 

C3  =  (P / a)[6 / (Ii  +  1  +  6)]  «  (26) 

Vo  a£Qi 

If  the  Type-II  biquad  is  used,  K  satisfies  the  following  equation 


K 4  +  (4  y  +  26)  K3  +  (6y2  +  6  yS  +  62)K2  +  (4  y3  +  6y26  +  2  yS2  +  —)K 


a 


,4  ,  <->„.3  c  ,  „,2  c2  ,  yPS  +  ^-SC 


Hv  +  2  yaS  +  y *6*  + 


0. 


a 


In  the  same  way,  it  can  be  seen  that 


K  «  y<5/(fi0T)  -  1 

C4  =  (^)/[Q(/f  +  l  +  «)(A'  +  l)]«2_. 

vo 

The  other  capacitor  values  can  now  be  computed  as 


(27) 


(28) 

(29) 


C5  =  k 

r  _  Csb2 

7(77  +  y  +  6)2{K  +  y) 
n  -  iC^K  +  y)(I\  +  y  +  6) 


(30) 

(31) 

(32) 
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Note  that  if  a  zero  outside  the  unit  circle  is  required,  then  one  more  capacitor  Cg  can  be  added  as 
shown  in  the  Fig.  3,  and  the  design  equations  can  also  be  obtained  in  the  same  way.  It  is  apparent 
from  the  above  analysis  that  the  use  of  the  charge  differencing  technique  can  reduce  the  required 
capacitor  spread  ratios  to  approximately  \J8  j(Q.0T).  This  means  that  a  comparable  saving  in  the 
silicon  areas  can  be  achieved  with  these  VLT  circuits.  For  instance,  in  a  lowpass  filter  design  (later 
used  in  the  implementation  of  the  parallel  filter  bank  with  Q o  =  0.707,  fl0  =  2iv  •  200 Hz,  and 
fs  =  125 kHz,  and  assume  unity  DC  gain,  y  =  1,  8  =  0.5),  the  capacitor  values  which  Type-II 
biquad  needed  are 

CA  =  CB  =  5.444  (33) 

C4  =  1-572  (34) 

C5  =  1.  (35) 

In  the  other  word,  the  total  capacitance  is  22.46  in  such  an  implementation.  If  compared  with  the 
Type-II  GOC  conventional  biquad  shown  in  Fig.  11  which  requires  the  total  capacitance  of  204.66, 
only  about  11%  of  capacitance  area  is  needed.  Even  compared  to  the  most  area-efficient  biquads  so 
far,  more  than  25%  capacitance  area  can  be  saved  easily[T2].  SWITCAP-II2  simulations  of  both  the 
ideal  case  (op-amp’s  gain=co)  and  the  non-ideal  case  (op-amp’s  gain=1000)  are  shown  in  Fig.  7. 
The  results  illustrate  the  robustness  of  these  area-efficient  VLT  biquads,  in  that  minimal  changes 
in  the  transfer  functions  occur  that  are  due  to  the  non-ideal  effects  of  the  op-amps. 

A  LPN  filter  has  also  been  fabricated  and  tested.  Both  the  frequency  and  time  responses  are 
shown  in  the  Figs.  8  and  9  respectively.  The  notch  frequency  deviation  from  its  designed  value 
is  less  than  1%  without  any  post-fabrication  tuning.  This  can  be  improved  further  if  larger  unit 
capacitor  is  used  instead  of  the  0.2pF  used  here.  A  single  tone  with  swept  frequency  was  used  to 
measure  the  entire  frequency  responses  of  the  fabricated-IC  chip.  The  responses  match  very  closely 
with  designs  already  discussed.  Furthermore,  the  DC  unity  gain  error  is  found  to  be  within  1%. 

3.3  Circuit  Design  II:  Conventional  GOC  Biquads 

As  we  mentioned  earlier,  conventional  switched-capacitor  circuits  have  a  capacitance  spread  of 
approximately  1/(D0 T).  This  ratio  is  not  problematic  in  the  high  frequency  channels.  Due  to  its 

2 SWITCAP-II  is  a  simulation  program  of  switched-capacitor  circuits. 
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Figure  7:  Simulation  result  of  the  area-efficient  biquad. 


Fiea  Rcsoomc  Ch  2  :  1  Avo=1G0 


Figure  8:  Measured  frequency  response  of  the  integrated  circuit  Type-II  area-efficient  L.PN  biquad. 
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X:5.50842mSec  V: -1.0254V 


X:4.776mSec  Y:S2243mV 


Figure  9:  Measured  time  domain  response  of  the  integrated  circuit  Type-II  area-efficient  LPN 
biquad  (input  frequency:  200 Hz,  input  signal:  lower  trace,  and  output  signal:  upper  trace). 


simpler  structure  than  VLT  biquads,  these  kind  of  circuits  are  preferred  in  the  applications  of  high 
frequency  cochlear  channels.  The  conventional  Type-I  and  Type-II  GOC  biquads  are  designed  and 


shown  in  in  Figs.  10  and  11,  respectively.  The  transfer  function  is  given  by 
H(s  Cs  ac  1  -  (2  -  7C  -  T]c)z~l  +  (1  -  t]c)z-2 

C2  7c  1  -  (2  -  ac  -  Pc)z~l  +  (1  -  Pc)z~2 

where, 


For  the  Type-I  biquad, 


For  the  Type-II  biquad. 


7c  = 
Vc  = 


C1C5 


cAc6 

c\c7 

cAc6 


Ci  C  2 
CaCb 
C!C3 
CACB ' 


CxC2 

CACg  +  CAC± 
C4C4 

CACB  +  C4C4 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 
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Figure  10:  Type-I  GOC  conventional  biquad 


Figure  11:  Type-II  GOC  conventional  biquad. 


Following  the  same  derivation  as  that  described  in  Section  3.2,  and  assigning  C\  =  Ci  -  C0 1  = 
C02  =  1,  we  obtain  the  following  capacitor  values  for  the  Type-I  biquad 


tn 

O 

II 

II 

sTr^i 

ii 

1 

%  T 

(43) 

a 

1  1 

(44) 

Qo  flo T 

In  the  same  way,  if  the  Type-II  biquad  is  used,  the  capacitor  values  are 

cA  =  Cb  = 

1 

Q.qT 

(45) 

C  -  P  ~ 
4  ~ 

1 

Qo 

(46) 

Where, 

Cs  = 

k 

(47) 

II 

O 

C5 

lCA 

(48) 

cv  = 

nCe,cA- 

(49) 

As  mentioned  earlier,  the  analysis  illustrates  that  Type-I  biquad  is  suitable  for  high  Q  realizations, 
while  Type- II  biquads  are  more  appropriate  for  low  Q  filters.  As  an  example,  using  a  Type-II  biquad 
to  implement  the  earlier  described  lowpass  filter  but  with  a  higher  pole  frequency  (fi0  =  27T-6.4 kHz), 
we  can  obtain  the  simulation  results  of  Fig.  12. 


3.4  Circuit  Design  III:  Highpass  Filter 


The  function  of  the  highpass  filter  is  to  pre-emphasize  the  input  signal,  and  hence  the  corner 
frequency  of  this  filter  should  be  high  enough  to  cover  the  entire  signal  bandwidth  of  interest.  A 
pole  frequency  of  flo  =  27r  •  12.8 kHz  is  used  here.  The  GOC  first  order  highpass  filter  which  employs 
the  offset  storage  capacitor  is  designed  and  shown  in  Fig.  13.  The  transfer  function  is  given  by 

{Ca  +  C2)-Caz-'- 


Considering  the  highpass  transfer  function  in  the  s-domain  where  k  is  the  gain  factor  and  fl0  is  the 
pole  frequency, 


H(s) 


~k 


s/ Qq 
1  +  s/£Iq 


(51) 
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Frequency  (Hz) 


Figure  12:  Simulation  result  of  the  conventional  GOC  biquad. 


Figure  14:  Simulation  result  of  the  conventional  GOC  first  order  highpass  filter. 


Performing  the  bilinear  transformation,  the  transfer  function  in  the  z-domain  is  obtained  as 

A  (1  -  z~l)x 0 

j  (1  +  x0)  —  (x0  —  1)2_1  ’ 

where,  Xq  —  2/(Q.oT),  and  T  is  the  sampling  period.  Comparing  the  coefficients  of  eqns.  50  and  52, 

and  choosing  the  capacitor  Ci  as  the  unit  capacitor,  we  obtain 


CA  =  (53) 

Cl  =  kj.  (54) 


Because  of  the  relative  insensitivity  of  the  circuit  responses,  C0  and  Cm  can  also  be  chosen  as  the 
unit  capacitors.  SWITCAP-II  simulations  of  the  pre-emphasis  filter  with  gain  k  =  5  are  in  Fig.  14. 
The  unity  gain  is  around  the  2.8  kHz  in  this  case. 


3.5  Circuit  Design  IV:  Sum-Gain  Amplifier 

In  a  parallel  filter  bank  in  which  biquads  elements  are  shared  among  neighboring  channels,  it 
is  essential  to  have  a  sum-gain  amplifier  which  can  sum  the  output  signals  of  several  biquads 
simultaneously.  A  new  sum-gain  amplifier  (SGA-SI)  which  permits  the  inputs  to  be  sampled  with 
the  same  phase  is  proposed  in  Fig.  15.  The  operational  principle  of  this  SGA-SI  can  be  stated  as 
follows:  Consider  input  V)„i  first,  while  V,'„2  and  Cm3  are  set  to  zero.  In  phase  2,  Ca 2  is  discharged. 
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Figure  15:  Sum-gain,  amplifier  with  same  input  phases  (SGA-SI). 


But,  since  in  phase  1,  its  left  plate  is  connected  to  the  inverting  terminal  of  the  op-amp  only,  no 
charge  can  be  injected  into  Ca 2,  and  the  output  remains  at  ground  level.  This  helps  Ca\  to  discharge 
in  this  phase.  During  phase  2,  C\,  Ca  1,  and  the  op-amp  constitute  a  non-inverting  gain  stage,  which 
gives  Vout(n  +  1/2)  =  (C1/C0i)Vjnl(n).  The  input  Vm3  can  be  analyzed  in  similar  fashion.  Next, 
set  Vini  and  Vj„3  to  zero.  In  phase  2,  Ca 2  is  discharged  as  in  the  previous  case.  But  in  phase  1,  C 2, 
C, 12,  and  the  op-amp  constitute  an  inverting  ampbfier.  Hence,  Vout(n)  =  —(C2/Ca2)Vin2(n).  This 
voltage  is  stored  in  Ca  1,  and  is  to  be  held  constant  during  the  next  phase  2.  The  final  output  is 
the  superposition  of  all  the  3  inputs.  The  transfer  function  is  given  by 


Vout(z) 


l^-Vinl (z)  -  ^-Vm2(z)  +  p-Vi,a 

^  al  ^  a2  ^  al 


Hence,  the  output  is  the  sum  of  Vjni,  V{n 3,  and  — Vjn 2,  all  scaled  by  appropriate  gain  factors. 
The  gain  factors  can  be  individually  controlled  through  C\,  C2,  and  C3.  Moreover,  this  sum-gain 
amplifier  can  be  designed  to  include  any  number  of  inputs  for  either  addition  or  subtraction.  One 
important  factor  which  needs  to  be  considered  carefully  is  the  slew  rate  because  the  op-amp  is  reset 
in  one  phase.  Taking  this  design  for  example,  the  clock  period  is  8 fis,  and  voltage  level  is  5.  The 
slew  rate  is  about  1.25  V/ ps  and  this  can  be  easily  achieved.  Even  so,  the  circuit  should  not  be 
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Figure  16:  Simulation  result  of  sum-gain  amplifier  with  same  input  phases  (SGA-SI). 

used  in  high  frequency  applications  because  of  the  low  effective  op-amp  gain.  Simulation  results 
with  input  V, n3  =  0  are  shown  in  Fig.  16.  In  addition,  The  responses  of  the  fabricated  sum-gain 
amplifiers  are  shown  in  Figs.  17  and  18.  The  testing  results  completely  satisfy  our  designs. 

3.6  Simulation  Results  of  the  Parallel  Filter  Bank 

Combining  all  individual  circuits  described  above,  we  have  implemented  the  shared  parallel  filter 
bank  shown  earlier  in  Fig.  1.  To  demonstrate  the  validity  of  these  designs,  the  parallel  filter  bank 
(including  the  pre-emphasis  highpass  filter)  was  simulated  by  SWITCAP-II  and  the  responses  are 
shown  for  4  channels  in  Fig.  19.  Note  that  the  channel  transfer  functions  are  strictly  translated 
relative  to  each  other  with  a  constant  Q.  The  change  in  the  relative  amplitudes  of  the  outputs  is 
due  to  the  pre-emphasis  filter.  This  architecture  saves  at  least  25%  of  the  silicon  area  needed  in 
the  non-sharing  topology. 
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Figure  17:  Measured  response  of  non-inverting  channel  of  the  integrated  circuit  SGA-SI  with  input 
frequency  =  IkHz  and  gain  =  9.475  (input  signal:  upper  trace,  and  output  signal:  lower  trace). 


Inst  Time  Chan  1 


Inst  Time  Chan  2 


X:1.8859mSec  Y:1. 7729V 

Figure  18:  Measured  response  of  inverting  channel  response  of  the  integrated  circuit  SGA-SI  with 
input  frequency  —  1  kHz  and  gain  =  18.8  (input  signal:  upper  trace,  and  output  signal:  lower 
trace). 
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1.2 


Figure  19:  Simulation  result  of  auditory  wavelet  transform  using  a  parallel  filter  bank  (including 
the  pre-emphasis  filter). 

4  Cascade  Filter  Bank 

4.1  System  Design 

The  other  possible  architecture  to  implement  wavelet  transforms  is  the  cascade  filter  bank.  The 
system  topology  of  a  cascade  filter  bank  implementation  is  shown  in  Fig.  20.  This  filter  bank  is 
made  up  of  a  first  order  pre-emphasis  highpass  filter  followed  by  several  cascaded  second  order 
lowpass  filters.  The  advantage  of  this  architecture  is  that  it  requires  a  smaller  silicon  area  when 
compared  to  the  parallel  architecture.  If  the  pole  frequencies  of  the  lowpass  filter  stages  are  chosen 
in  an  exponential  relationship,  the  lowpass  filters  become  dilated  with  respect  to  one  another.  But, 
the  dilation  relationship  among  the  channels  is  no  longer  exactly  preserved  when  the  lowpass  filters 
are  cascaded.  For  instance,  it  can  be  expected  that  the  last  channel  in  the  chain  should  have 
sharper  roll-offs  than  earlier  channels.  However,  this  may  not  be  a  problem  in  some  applications 
since  most  of  the  deviations  from  a  strict  wavelet  transform  occur  in  the  stop-bands  of  the  filter 
transfer  functions.  Finally,  another  important  feature  of  the  cascade  architecture  is  that  high  gains 
can  be  achieved  by  combining  modest  gains  of  many  stages  like  the  pseudoresonance[3].  Thus,  only 
small  Q  factors  are  needed  compared  to  the  values  typically  needed  in  the  parallel  filter  bank. 
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Figure  20:  System  topology  of  auditory  wavelet  transform  using  a  cascade  filter  bank. 

4.2  Circuit  Design 

The  circuits  designed  in  Section  3  can  also  lie  used  to  implement  the  cascade  filter  bank  by  modifying 
the  sign  of  the  transfer  function.  Since  only  lowpass  filters  are  used,  the  forward  path  capacitors 
Ce,  and  C 7  of  Figs.  3,  4,  10,  and  11  can  be  omitted.  To  obtain  positive  transfer  function,  we  take 
advantage  of  certain  SC  circuit  properties.  For  instance,  it  is  well  known  that  by  changing  the 
switching  phases  associated  with  a  capacitor,  either  positive  or  negative  resistor  can  be  simulated. 
Thus,  interchanging  the  switching  phases  of  LHS  plate  of  capacitor  C 5,  the  positive  transfer  function 
can  be  obtained  easily  without  requiring  an  additional  op-amp,  and  the  design  equations  are  the 
same  as  those  described  in  Section  3.  One  of  the  important  considerations  when  cascading  SC 
circuits  together  is  the  ripple  delay  since  it  may  make  it  necessary  to  use  high  speed  op-amps.  In 
our  biquad  circuit  design,  the  outputs  of  op-amps  are  sampled  in  the  different  phases,  hence  no 
ripple  delay  will  propagate  across  the  biquads  such  that  no  fast  op-amps  are  required.  By  using 
the  same  biquad  circuit  topology  with  a  different  clocking  scheme  (interchanging  switches  ‘1’  and 
‘2’)  as  the  adjacent  lowpass  filters,  one  can  cascade  these  filters  together  easily. 
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Figure  21:  Simulation  result  of  auditory  wavelet  transform  using  a  cascade  filter  bank  (including 
the  pre-emphasis  filter). 

4.3  Simulation  Results  of  the  Cascade  Filter  Bank 

To  demonstrate  the  validity  of  such  a  system  design,  and  in  order  to  compare  it  to  the  parallel 
architecture,  the  cascade  filter  bank  including  the  pre-emphasis  highpass  filter  was  simulated  by 
SWITCAP-II.  The  responses  are  shown  in  Fig.  21.  Note  that,  the  capacitance  spread  ratio  of  less 
than  10  is  needed  in  this  example  confirming  the  circuit  area-efficiency  which  we  mentioned  earlier. 
Since  no  sum-gain  amplifiers  are  required  here,  the  entire  silicon  area  is  only  about  64%  of  the 
parallel  filter  bank. 

5  Conclusions  and  Discussions 

We  have  discussed  in  this  paper  the  hardware  design  and  realization  of  auditory  wavelet  trans¬ 
forms  using  switched-capacitor  circuits.  The  desired  transfer  functions  could  either  be  systemat¬ 
ically  decomposed  into  rational  transfer  functions  by  the  recently-developed  wavelet  system  (WS) 
technique[13,  14]  as  shown  in  Appendix  A,  or  approximated  by  the  proposed  dilating-functions 
filter  bank  to  allow  for  circuit-sharing.  It  was  also  noted  in  Appendix  A  that  in  some  cases  the 
dilations  inherent  WS  approximations  may  also  be  exploited  in  the  design  of  shared-component 
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parallel  filter  banks.  In  this  paper,  lowpass  transfer  functions  were  chosen  as  the  desired  filter 
shapes  so  as  to  facilitate  the  comparison  between  two  possible  architectures  for  implementing  the 
wavelet  transforms,  namely  as  a  parallel  or  as  a  cascade  filter  bank.  1 

The  hardware  implementation  was  made  possible  by  several  new  circuit  designs.  Specifically, 
extremely  area-efficient  designs  were  developed  to  implement  very  large  time-constant  filters.  The 
designs  employ  a  new  charge  differencing  (CD)  technique  to  reduce  significantly  the  capacitance 
spread  ratios  needed  in  the  filter  banks.  Also,  a  new  sum-gain  amplifier  (SGA-SI)  was  designed 
which  permits  the  inputs  to  be  sampled  with  the  same  phase,  thus  facilitating  the  sharing  of  several 
biquads  among  neighboring  channels  in  the  parallel  filter  bank.  All  the  proposed  circuits  have  been 
successfully  fabricated  and  tested.  The  precise  responses  of  the  IC  chips  show  the  advantages  of  the 
new  designs.  Furthermore,  the  circuit  design  equations  were  fully  derived  in  order  that  arbitrary 
future  designs  of  the  wavelet  transform  filters  can  be  implemented. 

Based  on  these  circuits,  both  the  parallel  and  cascade  filter  banks  were  simulated  (Figs.  19  and 
21).  The  simulation  illustrated  that  both  filter  banks  have  similar  robust  response  characteristics. 
One  difference  between  these  two  is  that  faster  decaying  slopes  seen  in  the  stop-bands  of  the  cascade 
realization.  Since  the  signal  is  very  small  in  the  stop  band,  this  difference  is  generally  not  critical 
in  computing  the  wavelet  transform.  The  cascade  filters  also  needed  smaller  silicon  areas  and  lower 
Q  circuits  when  compared  to  the  parallel  realization.  However,  more  robust  circuits  with  low  DC 
offset  voltages  must  be  used  in  the  cascade  realization. 

The  other  key  feature  of  SC  circuits  from  the  viewpoint  of  wavelet  transforms  is  that  by  changing 
the  sampling  frequency  by  a  factor  of  u,  all  zero  and  pole  frequencies  can  be  changed  by  the 
same  factor.  Hence,  by  connecting  several  similar  chips  together,  each  with  a  different  sampling 
frequency,  more  filter  channels  can  be  readily  implemented.  Besides,  the  dilation  constant  a  can 
be  very  precisely  controlled  by  either  the  clock  frequency  or  the  capacitor  ratios  in  SC  circuits.  For 
instance,  in  the  many  applications  of  the  wavelet  transforms  where  the  dilation  constant  a  is  equal 
to  2,  simple  frequency  division  digital  circuits  can  be  used. 
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Appendix 


A  Rational  Wavelet  System  Approximations 

Rational  WS  approximations  to  transfer  functions  are  based  on  wavelet  decompositions  of  the 
Hardy  space  H2(n+),  where  n+is  the  open  right-half  complex  plane  5ie  s  >  0,  in  the  case  where 
the  analyzing  wavelet  is  chosen  to  be  rational. 

Definition  A.l  A  function  F  which  is  analytic  in  n+zs  said  to  belong  to  the  class  H2(n+)  if, 

sup  [  \F(x  +  iy) j2  dy  <  oo.  '  (56) 

x>0 

By  the  Paley-Wiener  theorem,  H2(n+)  is  the  image  of  L2(0,co)  under  the  Laplace  transform. 
Hence  H2(n+)  contains  transfer  functions  of  causal,  linear,  time-invariant  systems  with  square- 
integrable  impulse  responses.  We  adopt  the  following  notation. 

RH2(n+)  :  Real-rational  functions  in  H2(n+). 

H^(n+)  :  Functions  in  II2(n+)  which  are  Laplace  transforms  of  real-valued  functions  in  L2(0,  oo)  .. 

WS  approximations  arise  from  specific  truncations  of  decompositions  of  H2(n+),  via  dilations 
and  complex  translations  of  a  single  real-rational  analyzing  wavelet.  The  set  of  dilations  and 
translations  of  the  real-rational  analyzing  wavelet  are  chosen  so  as  to  form  a  frame  for  H2(n+). 
Frames  are  generalizations  of  orthonormal  bases,  which  are  defined  below.  We  refer  to  [20]  for 
further  details. 

Definition  A. 2  Given  a  Hilbei't  space  hi,  a  sequence  of  vectors  {hn}n  C  ' H ,  is  called  a  frame  if 
there  exist  constants  A  >  0  and  B  <  oo  such  that 

A||/||2  <£|</,/i„>  |2  <5||/||2,  (57) 

n 

for  every  f  £  hi.  A  and  B  are  called  the  frame  bounds. 

Given  a  frame  {hn}n  for  a  Hilbert  space  hi,  any  /  6  ’H,  may  be  represented  as 

/  =  £(/’  K)  S~lhn  =  £  (/,  s-X)  hn , 

n  n 

where  the  frame  operator  S,  is  defined  by  Sf  =  (/,  hTl)  hn.  The  following  theorem  summarizes 

the  main  ideas  behind  WS  approximations. 
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Theorem  A.l  ([14])  Let  $  £  RH2(II+)  be  an  admissible  analyzing  wavelet,  and  let  Oq  >  0,  bo 
be  such  that  (^,a0,b0)  generates  an  affine  frame  for  H2(II+),  i.e.  {^miTl}  =  '(a'ff  ■  -m60)j 

is  a  frame  for  H2(II+).  Let  S  be  the  associate  frame  operator.  Then,  any  F  in  H^IT*")  may  be 
represented  as, 

OO 

F  =  EE  (58) 

m  n= 0 

where,  each  Fm,n  (£  RH2(II+))  is  defined  by, 

Fm,n  _  ^F,S~  1«m,n)  ^m,n  +  (F,  S"1  *m,n) ¥m,_n ,  TTl  £  Z,  U  €  Z+  \{0} 

Fm,o  =  ^>lSri$mi0^mt0  meZ.  ■  (59) 


The  right  hand  side  of  eqn.  58,  is  referred  to  as  a  wavelet  system  (WS)  decomposition  of  F  £  H2(II+). 
WS  decompositions  represent  systems  with  transfer  functions  in  H2(II+),  via  infinite  sums  of  time- 
frequency  localized  finite-dimensional  systems.  Finite  truncations  of  WS  decompositions  give  real- 
rational  approximations  to  transfer  functions  in  H2(IT+).  WS  techniques  are  particularly  well-suited 
to  rational  approximation  of  transfer  functions  which  are  well-localized  in  time-frequency  (as  is  the 
case  for  the  cochlear  filters).  Such  finite  truncations  are  also  useful  as  ‘linear-in-parameters’  black¬ 
box  models  for  system  identification  (c.f.  [14]). 


A.l  WS  Approximation  of  Cochlear  Filters 


For  approximation  of  the  cochlear  filter,  we  use  the  following  H2(II+)  analyzing  wavelet  HI; 


1 

(s  +  7)2  -K2’ 


>  0, 


with  7  =  5.0,  £  =  1.0.  It  is  easily  verified  that  'i(s)  is  in  H2(II+)  and  furthermore  $  is  an  admissible 
analyzing  wavelet  i.e.  /^^(x  +  iy)dy  =  0  for  x  >  0.  It  can  also  be  shown  (see  [14])  that  for 
a0  =  2,  0  <  60  <  16.5,  (4/,  a0,  fro),  generates  an  affine  frame  for  H2(II+). 

The  specific  cochlear  filter  shape  of  interest  here  (see  Figure  23)  takes  the  form  of  an  unparam¬ 
eterized  model  derived  from  a  combination  of  empirical  data  and  mathematical  modeling  of  the 
basilar  membrane  [21].  Figure  22  is  a  plot  of  the  magnitudes  of  the  wavelet  expansion  coefficients 
{(F,  5,-1'Pm,n)}  (using  a0  =  2,  bo  =  1.0375).  A  key  feature  of  the  decomposition  shown  in  Figure 
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Figure  22:  Magnitude  of  wavelet  expansion  coefficients  for  the  cochlear  filter  transfer  function. 

22  is  the  large  number  of  zero  or  negligible  coefficients.  This  compactness  of  representation,  which 
arises  due  to  the  time-frequency  localization  of  both  the  cochlear  filter  transfer  function  and  the  an¬ 
alyzing  functions  (wavelets),  allows  us  to  construct  low-order  rational  approximants  to  the  cochlear 
filters.  Figure  23  shows  the  frequency  response  magnitude  of  a  12th  order  rational  approximation  to 
the  cochlear  filter  response  in  which  three  terms  ({(m, «)}  =  {(0,  7),  (1, 14),  (2,  28)})  from  the  WS 
decomposition  are  used.  Note  that  the  sharp  cutoff  on  the  high-frequency  side  is  well-approximated 
by  this  scheme. 

Remark:  In  general  WS  decompositions  serve  to  provide  rational  approximations  which  are  useful 
for  independent  implementation  of  filters  in  a  parallel  filter  bank.  However,  the  fact  that  the 
individual  terms  of  a  WS  decomposition  involve  dilations  of  the  analyzing  wavelet  may  in  some 
cases  also  be  exploited  in  designing  the  economical  shared-component  form  of  a  parallel  filter  bank 
implementation  described  in  Section  3. 
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Figure  23:  Frequncy  response  magnitude  of:  (i)  unparameterized  cochlear  filter  model  (dashed 
line),  (ii)  12th-order  rational  WS  approximation  to  cochlear  filter  (solid  line). 
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